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Chiral dynamics makes definitive predictions for the electromagnetic polarizabilities of hadrons 
near the chiral limit; but, agreement with experiment is tenuous in some cases. We provide an 
overview of lattice QCD methods to compute the electric and magnetic polarizabilities of hadrons. 
Central to these methods is the lattice simulation of quarks in uniform, classical electromagnetic 
fields. A long-term goal is the determination of polarizabilities directly from lattice computations, 
however, in the near term, one may need to rely on partially quenched chiral perturbation theory. 
Nonetheless the same striking predictions for the pion mass dependence of electric and magnetic 
polarizabilities can be made from chiral dynamics, and tested with lattice QCD. A particular focus 
is a novel new method to handle charged particle correlation functions in magnetic fields. 
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1. Electromagnetic Polarizabilities 

The electric and magnetic polarizabilities provide an opportunity for a stringent test of the 
chiral dynamics inside hadrons. Essential to the stoiy of chiral dynamics is the spontaneous break- 
down of chiral symmetry. For the case of two massless quark flavors, the QCD functional integral 
possesses an SU (2)l®SU (2)r symmetry, which is broken to the vector subgroup by the vacuum 
expectation value of the chiral condensate, < Y l Yr 0. In this picture, the pions emerge as 
Goldstone modes, and have masses that vanish in the absence of explicit chiral symmetry breaking, 
that is m\ = 0+m q \ < > \F 2 , where the explicit chiral symmetry breaking is parameterized by 
m q , the quark mass. This picture is an effective description of low-energy QCD provided the quark 
masses are small compared to the scale of strong interactions. In terms of hadronic parameters 
that are free from QCD renormalization scale and scheme dependence, this condition translates 
into m n /{A%F) <C 1; and, when met, allows for the quark mass dependence of low-energy QCD 
observables to determined systematically in an expansion about the chiral limit, m q = 0. 

In chiral perturbation theory, the interactions of hadrons with pions are constrained by the 
form of spontaneous and explicit chiral symmetry breaking. As a result, the bare hadron fields are 
dressed with pions in a Fock state expansion, schematically of the form 

|n°) = c \x°)+c l \x°x + x-) + ..., \N) = c \n)+ c:\px~) + ..., (1.1) 

where the bare fields are denoted with lower-case lettering, co's are related to the wave-function 
renormalization, and ci's are determined by the interactions of the theory. Depending on kinemat- 
ics, the higher Fock components may only be virtual states. Above we have chosen states that are 
electrically neutral to emphasize that dynamical fluctuations produce Fock components containing 
charged hadrons. The electromagnetic interaction thus serves as a probe sensitive to the higher 
Fock components of the bare hadron fields. Within chiral perturbation theory, the electric and mag- 
netic polarizabilities (cfe and j3 M , respectively) are determined to leading order entirely from higher 
Fock components. For the poins |jl|], and nucleons there is a characteristic singularity in the 
chiral limit of these quantities, 

;U (1-2) 



Experimental determination of polarizabilities can be achieved through the analysis of low-energy 
Compton scattering data. In principle, this is easiest for the proton, however, the low-energy limit 
is dominated by the total charge interaction, the Thomson cross section; and, while increasing the 
energy leads to increased sensitivity to polarizabilities, it also introduces higher-order response of 
the nucleon. For the most recent comprehensive analysis of proton Compton scattering, see [|3|]. 
For the neutron, one must use Compton scattering off deuterium or other light nuclei to extract 
polarizabilities, which introduces additional theoretical uncertainty. For pions, experiments have 
resorted to photo-pion production off the nucleon [Q], and new results are anticipated from recent 
COMPASS measurements using pion scattering off Primakoff photons |@j. 

2. Methods for Particles in Electric Fields 

The electromagnetic polarizabilities provide an opportunity for lattice QCD computations, as 
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they are quantities that have been subject to debate. A lattice determination of the polarizabilities of 
the deuteron, for example, would be a major contribution, as one would be able to study few-body 
dynamics in addition to the chiral behavior. Needless to say, lattice QCD is not yet at the point 
of such studies, however, there has been progress in treating strong interactions in the presence 
of uniform, classical electromagnetic fields, see [Q] for overviews. One might think the natural 
starting point for lattice computations of polarizabilities would be the Compton scattering tensor, 

T^ v (k',k) = J d\d 4 ye- ik! ^ +ik ^(H\T{J^ X )J v (y)}\H), (2.1) 



however, this presents two major complications. The low-energy limit in Eq. 2.1 is constrained by 
spatial momentum quantization conditions resulting from periodic boundary conditions imposed 
on the quark fields. Progress has been made on this front. The desired terms at second order in the 
photon frequency can be isolated by studying zero momentum derivatives of the quark propagators, 
and these can be computed approximately by using partially twisted boundary conditions, and 
taking the limit of vanishing twist angle This procedure leads one to the computation of certain 
hadronic four-point functions for which nearby intermediate states between the current insertions 
present an essential complication. While there has been progress in computing four-point functions 
in the case of the mesons, namely the Kl-K$ mass difference [g], it is likely that such methods will 
not be practicable for nucleons. 

The external field approach provides an alternate method to access electromagnetic polariz- 
abilities. One adds a classical electromagnetic field to QCD computations, and studies the sub- 
sequent external field dependence of hadronic correlation functions. Such dependence gives one 
access to hadronic couplings to the external field. To include an external electromagnetic field, one 
appends U (l)-valued links to the Si/ (3) color gauge links 



U^x)^U^ m ( X )U^( X ). (2.2) 

This must be done for valence quark propagators, and in the quark determinant used to generate 
gauge ensembles. The latter encompasses contributions due to the electric charges of sea quarks, 
and has only been achieved on lattices studying thermodynamics with staggered quarks. In the 
near term, calculations in weak external fields will exclude contributions from the sea quarks in 
gauge field generation. Gauge ensemble re-weighting techniques are a promising way to include 
sea quark charges [Qj. The light quark mass regime is problematic for the quenching of sea quark 
charges altogether, due to the so-called exceptionality of gauge configurations. Exceptional config- 
urations create an essential roadblock for post-multiplying electromagnetic links to existing gluon 
gauge configurations. Inasmuch as such configurations are not encountered, one can address the 
quenching of sea quark electric charges using chiral perturbation theory, for a discussion see 
and predictions exist for pion and nucleon polarizabilities as a function of the sea quark electric 
charges [ffly. As our methods are addressed with exploratory lattice studies, we couple the external 
fields to valence quarks only. 

A final word on the inclusion of uniform, external electromagnetic fields on a lattice. The peri- 
odicity of lattice quark fields leads to a quantization condition on the strength of external fields [|j~2|]. 
because the hyper-torus forms a closed surface that does not leak any flux. For electric and mag- 
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netic fields, the quantization conditions are of the form 1 

qg = 2xn/LP, qB = 27ln/L 2 , (2.3) 

where L is the length of a spatial direction, with all three spatial directions assumed to be of the 
same length, and j3 is the length of the temporal direction of the lattice. We write £ for the 
electric field to make clear that we are in Euclidean space, where the action density would appear 
as \F^ V F^ V = \ {b 1 + $ . Analytic continuation is necessary for results in Minkowski space; but, 
as we are interested in quantities that are perturbative in the strength of the field, the continuation 
is trivial. The Schwinger mechanism [JT3J] , which is a non-perturbative phenomenon, is fortunately 
absent in Euclidean space. In this section, we concern ourselves with external electric fields. 

2.1 Neutral Particles 

For neutral hadrons, the method is simple to explain. One includes an external electric field, 
and measures two-point correlation functions of hadrons. In the long Euclidean time limit, these 
correlation functions should follow an exponential falloff, 

Gg{%)=Zge- E ^\ (2.4) 

where E{$) is the energy of the hadron in the external electric field. For a neutron in an external 
electric field, the energy is given by [ |l6| ] 



E{S) = M n + (a E - At,?/4M„ 2 ) g 2 /2 + ..., (2.5) 

where terms of order <# 4 have been dropped. The contribution involving the square of the neutron 
magnetic moment arises from treating the neutron spin relativistically. In terms of neutron Compton 
scattering, the analogous contribution appears as a Born term, namely the second-order term arising 
from two interactions of leading-order. The exponential falloff of neutron correlators leads to an 
extraction of the quantity in parentheses. If one is interested in isolating the electric polarizability, 
one requires a method to determine the neutron magnetic moment, and this can be achieved by 



looking at the amplitudes in off-diagonal spin components of the correlator [|16|. In that study, an 



ensemble of anisotropic clover fermions [17] was successfully used to demonstrate the technique. 



The analogue of the Born term was shown to affect extraction of the electric polarizability by 50%. 
2.2 Charged Particles 

The study of charged particles in external electric fields is obviously complicated, however, 
the same philosophy can be applied. One determines hadronic correlation functions for various 
field strengths, and then matches onto the behavior expected from a single-particle effective action. 
This behavior is not a simple exponential falloff in Euclidean time. For example, the charged pion 
propagator should behave as 



G,(T) =Z g f * ■ e -W***Q**-\n^ (2 . 6) 
Jo v suing©* 



1 Strictly speaking such quantization conditions do not lead to just uniform electromagnetic fields. There are also 
gauge-invariant, finite volume artifacts involving the non-trivial holonomy of the external field JTit. While such contri- 
butions should be exponentially suppressed, ~ e^"'" L , lattice results suggest that the effect might be non-negligible even 
at large pion masses KM- That study employs methods different than those outlined here. 
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with E n {S) = m n + \(Xe$ 2 H ■ While fits to the non-standard T-behavior can be challenging, 

the technique has been successfully demonstrated in an exploratory lattice calculation of charged 



pion and kaon correlation functions Q19|]. Similar success was also achieved in the study of proton 



correlators using a generalization of the method to spin-half particles [|16|]. 

3. A Method for Charged Particles in Magnetic Fields 

The quantization condition for uniform magnetic fields has so far proven restrictive for the 
study of perturbatively small effects for hadrons. With increased lattice volumes, smaller magnetic 
field strengths can be accessed. The energy eigenstates of charged particles in external magnetic 
fields are described by Landau levels. For the charged pion, the long-time limit of the standard 
lattice correlation function should produce the standard exponential behavior 

G B (T)=£(7T(x,T) 7 r t (0,0)) =Z B e- E ^ + --, (3.1) 

x 

where Eq is the pion energy in the lowest Landau level, Eq(B) = m K + i^j- + \$mE> 2 H ■ The 

omitted terms in the long-time limit of the correlation function include excited hadronic states, as 
well as the higher Landau levels. For large values of the external magnetic field, the Landau levels 
will be widely separated in energy, and only the lowest Landau level will survive the long-time 
limit. As we are interested in perturbatively small magnetic fields, however, the narrow Landau 
level spacing, AE/M = \ QB\/M 2 , will lead to a pileup. For smaller values of the magnetic field, 
one will require longer times to separate out the contribution from the lowest Landau level. 

This complication can be sidestepped altogether [20]. Inspecting Eq. [3.1| , we see that the sum 



over all lattice sites projects the correlator onto zero spatial momentum, p = 0. In the presence of 
a magnetic field, it is impossible for all components of three-momentum to remain good quantum 
numbers. Consequently the correlator contains all Landau levels. A more judicious choice of 
correlation function is given by 

%(T) = £v5W(w(x > T)wt(0,0)) ) (3.2) 

X 

where Xj/o(x) is the coordinate wave-function of the lowest Landau level. The long-time limit of 
the correlation function has the same exponential falloff, however, the first omitted terms are from 
higher lying hadronic states just as in the absence of the magnetic field. This can be demonstrated 
with the Schwinger proper-time trick. Explicit projection of the lowest Landau level should be an 
economic technique in perturbatively small magnetic fields. 

As the technique must be practicable on the lattice, we investigated the effects of discretiza- 
tion on the lowest Landau level, and the effects of finite volume. For the former, we found that 



discretization effects on the wave-function of the lowest Landau level required in Eq. [^2] were neg- 
ligible. The discretization corrections to the energy of the lowest Landau level could have a more 
substantial effect competitive with the magnetic polarizability, 

E 2 (B,a) = ml + \QB\ - I -a 2 Q 2 +m x Pu) B 2 . (3.3) 



Finite volume corrections where shown to be important; but, can be treated as an application of the 
magnetic periodicity of the action [pi)]. 
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